Structural trends in clusters of quadrupolar spheres 
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The influence of quadrupolar interactions on the structure of small clusters is investigated by 
adding a point quadrupole of variable strength to the Lennard- Jones potential. Competition arises 
between sheet-like arrangements of the particles, favoured by the quadrupoles, and compact struc- 
tures, favoured by the isotropic Lennard- Jones attraction. Putative global potential energy minima 
are obtained for clusters of up to 25 particles using the basin-hopping algorithm. A number of struc- 
tural motifs and growth sequences emerge, including star-like structures, tubes, shells and sheets. 
The results are discussed in the context of colloidal self-assembly. 



I. INTRODUCTION 

The structure adopted by a collection of particles is ul- 
timately governed by the energetic interactions between 
the particles. It is therefore natural to ask what sort 
of structures are favoured by a given set of interactions, 
and also whether interactions can be chosen or manipu- 
lated in order to produce a particular target structure. 
The scope of both these questions is growing increasingly 
broad as it becomes possible to exert ever greater con- 
trol over the shape and form of the interactions between 
molecular and colloidal building blocks [1]. The motiva- 
tion for seeking deeper understanding in these areas is 
the desire to design novel materials and supramolecular 
structures with unusual and useful properties. 

Considerable structural variety is possible even for 
spherical particles and isotropic interparticle potentials. 
Simple van der Waals interactions between inert gas 
atoms promote near-spherical, highly-coordinated struc- 
tures, favouring icosahedral packing for clusters Q 
and close-packed crystals in the bulk. However, these 
structures can be suppressed by introducing a repulsive 
barrier into the pair potential at a distance close to y/2 
times the nearest neighbour separation [4]. In clusters, 
potentials of this form further promote icosahedral local 
order, but lead to polytetrahedral structures [5], or to 
less compact shapes that are either elongated or contain 
holes 0|. Local maxima in the pair potential can arise, 
for example, from the combination of short-range deple- 
tion attraction between colloidal particles with partially 
screened long-range Coulomb repulsion [7[ , in which case 
the accumulated charge of a growing cluster can effec- 
tively limit the size of aggregate that forms. It is also 
possible to favour less highly-coordinated order by careful 
design of isotropic pair potentials. For example, inverse 
statistical mechanical techniques can be used to derive 
isotropic potentials that favour crystalline but non-close- 
packed bulk structures, such as the diamond and wurtzite 
lattices 0|. 

A vast range of structures become possible when ei- 
ther non-spherical particles or anisotropic interactions 
are considered. Evolution has selected molecules with 
interactions that lead to the self-assembled structures 
that we observe in living matter, including pseudo 



one-dimensional filaments, two-dimensional membranes, 
tube-like channels and pores, and shell-like capsules. Fas- 
cinating examples in the latter category are the capsids 
of viruses, many of which are spheroidal and are con- 
structed from a specific number of copies of a small num- 
ber of different proteins [9[ . It has long been known that 
the capsomers of some viruses can assemble in vitro into 
empty shells even in the absence of the viral genetic ma- 
terial [TO], providing a powerful demonstration of how 
molecular shape and interactions dictate the structure of 
aggregates. For proteins, the symmetry and binding in 
oligomers is determined by the contacts between neigh- 
bours in the complex, and it is now becoming possible 
to engineer the quaternary structure by modifying the 
contact surface through mutations of the amino acid se- 
quence \u\ . 

Drawing inspiration from Nature, and encouraged by 
rapid advances in the synthesis of tailor-made building 
blocks, computational scientists are trying to understand 
the principles of self-assembly and how these can be used 
to build designed structures. Explicit models of polyhe- 
dral shell assembly have shown how difficult this objec- 
tive can be [l2[. While entropic and kinetic considera- 
tions are undoubtedly crucial for successful self-assembly, 
there is also a clear requirement for target structures to 
be energetically stable and kinetically accessible pjj 03 • 
A logical starting point is therefore to design building 
blocks with attractive sites or "patches" in geometries 
compatible with the target structure. Monodisperse dis- 
crete objects, as well as continuous structures, such as 
sheets, can be constructed in this way [TEL 

Highly directional interactions can also be achieved 
through a non-uniform charge distribution in particles 
that remain neutral overall. For example, dipolar par- 
ticles will endeavour to form chains, since the head-to- 
tail arrangement of two dipoles is low in energy and me- 
chanically stable. The tendency to form chains is par- 
tially frustrated in the presence of competition for more 
compact arrangements arising from an isotropic van der 
Waals or depletion attraction [l?], EH • A simple model 
incorporating these features is the Stockmayer potential, 
consisting of particles with a Lennard- Jones (LJ) site plus 
a central point dipole. By adjusting the relative strength 
of the LJ and dipolar contributions, the energetically 
most stable morphology of the 13-particle Stockmayer 
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cluster changes in four stages from a distorted icosahe- 
dron to a closed ring [I?], [lij]. For slightly larger sizes, 
knots, links and coils emerge [20( | . 

Van Workum and Douglas have investigated the self- 
assembly of chains in low-density Stockmayer fluids 
(Hf, regarding the process as a form of reversible 
polymerisation. These authors have also considered 
a natural extension of the Stockmayer potential to 
higher multipoles, in particular an LJ site plus point 
quadrupole [2l], [22[. Quadrupole-quadrupole interac- 
tions favour the formation of extended two-dimensional 
sheets, which can produce tubes when the edges become 
connected. The LJ-plus-multipole class of potentials 
therefore provides control over the preference for compact 
(three-dimensional), sheet-like (two-dimensional) and 
chain-based (one-dimensional) structure in self-assembly. 
These tendencies compete with each other. In the present 
contribution we provide a systematic survey of the struc- 
ture of small clusters of quadrupolar spheres in an at- 
tempt to identify and understand the structural motifs 
that emerge from the frustration between the isotropic 
and directional components of the potential. 



II. METHODS 

A. Model potential 

The quadrupolar sphere is modelled as an isotropic 
Lennard- Jones site with a point quadrupole of variable 
strength superimposed [U, [23[. The pair potential, 
which we denote as LJQ, is of the form 
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where R^ is the vector from particle i to j, Rij is the 
magnitude of this vector, and Qi represents the orienta- 
tional degrees of freedom of particle i. The parameters 
u and <j are the Lennard- Jones dimer equilibrium well 
depth and separation, respectively, and will be used as 
the units of energy and length henceforth. Vq is the 
quadrupole-quadrupole interaction, which depends on 
the component (s) of the quadrupole tensor involved. 

In this work, we consider the two quadrupolar arrange- 
ments of charges depicted in Fig. dJ In each case, the 
point quadrupole is reached by taking the limit in which 
the separation of the charges d goes to zero, while the 
strength Q = qd 2 of the quadrupole moment is held fixed. 
For the linear arrangement of charges, the interaction be- 
tween two point quadrupoles i and j can be written in 
terms of the unit vectors e iz and ej z along the body-fixed 
z-axes of the particles as 
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FIG. 1: Definition of the charge distributions in the "linear" 
(left) and "square" (right) quadrupoles with the local axis 
frame. The point quadrupole is obtained in each case by 
taking d — > while keeping Q — qd 2 fixed. 
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ej a • Hij/Rij (note the sign convention with respect to 
the direction of R^), with a and /3 representing y or 
z. The dimensionless parameter Q* = Q/(47re<j 5 ^) 1 / 2 is 
the reduced quadrupole strength, e being the dielectric 
permittivity of the medium. For the "linear" quadrupole 
defined by Eq.([2|), we represent the orientation of the 
quadrupole using spherical polar angles. 

In the d — » limit, the pair interaction for the square 
arrangement of charges in Fig. [1] is 
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Three variables are now required to specify the orienta- 
tion of the quadrupole, and we have chosen to represent 
the vectors e iy and e iz in terms of Euler angles. The 
LJQ potential with Vq = has received less atten- 

tion [2l|,|22[ in the past than LJQ with Vq = V^ in . 

The linear quadrupole corresponds directly to the 
spherical tensor component Q2O' The square quadrupole 
can be reached continuously from the linear arrange- 
ment by increasing the angle in Fig. [T]from to 7r/2, 
which corresponds to introducing a contribution from 
the component Q22C At = 7r/2, the combination is 
f Q20 + \ V3Q22c The interactions between components 
of the quadrupole in the spherical tensor representation 
are tabulated in Appendix F of Ref. [24| . 

We will need to perform local geometry optimisations 
using the LJQ potentials, and have therefore derived and 
coded their analytic derivatives with respect to the Carte- 
sian position coordinates and the angular orientational 
variables. 



B. Global optimisation 
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We performed unbiased searches for the global minima 
of clusters bound by the LJQ potentials using the basin- 
hopping algorithm [25[, in which a Monte Carlo simu- 
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lation is run on a transformed potential energy surface 
(PES) by performing a local minimisation of the energy 
at each step. The local minimisation [26] is key to the 
success of basin- hopping [27| , and is a feature shared by 
other efficient methods of global optimisation for clusters, 
such as certain genetic algorithms [281 ] . 

For a given number of particles, iV, and quadrupole 
strength, Q*, several runs seeded from different random 
initial positions and orientations were performed. The 
number of Monte Carlo steps required to find a putative 
global minimum reliably in independent runs depends 
strongly on the size of the cluster and the strength of 
the quadrupole moment. It is also important to select a 
reasonable temperature for the accept /reject step in the 
basin-hopping runs. Although the success of the method 
is not very sensitive to the temperature, it must be high 
enough for the search to escape from local traps, but not 
so high that we fail to sample the low-lying minima in 
each region sufficiently. A fixed reduced temperature of 
kT/u = 1 was often found to work well. 

To generate a structural map for the clusters, it is nec- 
essary to explore the two-dimensional parameter space 
defined by the size of the cluster and the strength of the 
quadrupole. The straightforward approach of running 
basin-hopping on a grid of Q* points for each N would be 
inefficient, since a small change in Q* will often lead only 
to a relaxation of the global minimum, with qualitative 
changes to a new structure occurring at larger intervals 
in Q*. We have therefore devised a surveying scheme 
with an iterative element, designed to identify the values 
of Q* where the identity of the global minimum changes 
for a given N. 

The algorithm begins with thorough searches for the 
global potential energy minimum at two values of the 
quadrupole strength, Q* ow and Q^igh' that are ^ ar en °ugh 
apart to lead to qualitatively different structures. These 
structures are then relaxed by local minimisation on a 
grid of Q* points that lie between Q* ow and Qh igh , re ~ 
suiting in the correlation of the energy of each structure 
with Q*. During this process, it is possible that one or 
both of the minima will disappear at some value of Q* 
due to a catastrophe in the PES [29j], leading to a sud- 
den change in energy as the structure falls into a different 
basin of attraction. In this case, a new basin-hopping run 
is performed at the Q* where the catastrophe occurred 
and the scan in Q* is continued. 

Eventually, the energies of the relaxed structures ini- 
tiated from Q* ow and Q£i g h cross a ^ some quadrupole 
strength QJross- New basin-hopping runs are then per- 
formed to identify the global minima at QJross =^ a 
little above and below the crossing point. If the basin- 
hopping run at Q cross — ^Q* returns the same structure 
and energy as the relaxed structure from Q* ow , then we 
assume that this structure was the global minimum not 
only at Q* ow and at Ql TOSS — 5Q*, but also at all values 
of Q* in the intervening range. In other words, we as- 
sume that there are no reentrant global minimum struc- 
tures. This assumption is not only intuitively reasonable, 



given that Q* continuously changes the potential from 
isotropic van der Waals attraction to highly directional 
electrostatic interactions, but it is also borne out by care- 
ful checks of particular cases. The range Q* ross + SQ* to 
Qhigh was treated analogously. Since a local relaxation is 
much faster than a full basin-hopping run, this procedure 
is far more efficient than using basin-hopping afresh at 
each intermediate Q* value. 

If, on the other hand, the basin-hopping runs at 
Qcross =^ $Q* return new structures with lower energy 
than the relaxed structures from Q* ow and Q£ igh , then 
this global minimum supersedes them and was in turn re- 
laxed at values of Q* in both directions away from Q* ross 
until the energy rose above those of the previous struc- 
tures. A new check for the true global minimum must 
now be performed at this crossing point. The procedure 
was terminated when no lower minima were found at the 
crossing points of relaxed structures. Hence, full basin- 
hopping runs need only be performed close to the loca- 
tions where the identity of the global minimum changes. 



III. RESULTS 
A. Local coordination of quadrupoles 

For a given separation R, the energetically optimal 
arrangement of two point quadrupoles is with the local 
y axis on one particle and the local z axis on the other 
aligned with R and the other two local axes coplanar. 
This is true for any value of in Fig. [TJ but we will refer 
to the arrangement as a "T-shape," which is most clearly 
seen for the linear case, = 0. 

For the trimers, denoted LQ3 and SQ3 for the linear 
and square quadrupoles, respectively, triangular arrange- 
ments are optimal, as shown in the upper panels of Fig. [21 
Despite the significant distortion away from three ideal 
T-shaped pair interactions, the energies of the trimers 
are each about 2.8 times the respective dimer energies at 
Q* = 5. 

In the tetramers (lower panels of Fig. [2]), the strain is 
relieved, making four undistorted T-shapes possible. For 
an interior angle of 45°, the slipped-parallel arrangement 
of the next-nearest neighbours is also favourable, further 
lowering the total energy of the tetramers to about 4.5 
times that of the respective dimers at Q* = 5. An impor- 
tant difference between the linear and square quadrupoles 
is demonstrated by the tetramers. The axial arrangement 
of charges in the linear quadrupole means that rotation 
of a quadrupole about a local z axis makes no differ- 
ence to the energy. LQ4 is therefore able to lower its 
energy by twisting the quadrupolar axes slightly. This 
distortion places diagonally opposite particles above and 
below the plane of the projection in the lower-left panel 
of Fig. O allowing next-nearest neighbours to approach 
more closely. In contrast, SQ4 does not have this flexi- 
bility, and the structure in the lower-right panel of Fig. [2] 
is completely planar. 
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FIG. 2: Characteristic coordination motifs for the linear (left) 
and square (right) quadrupoles: global minima for the trimer 
(upper panels) and tetramer (lower) for a strong quadrupole 
(large Q*). Extended charge arrangements are shown for il- 
lustration only; all calculations are in the point quadrupole 
limit. In the case of LQ4 (lower left), the quadrupole axes are 
not coplanar. 



We note that the T-shape and slipped-parallel pair ge- 
ometries are stationary points for dimers of the linear 
quadrupolar molecule carbon dioxide. However, in con- 
trast to the LJQ model, the T-shape of (002)2 is a saddle 
point, while the slipped-parallel geometry is stable [30[. 



B. Strong and weak quadrupole limits 

The T-shaped nearest-neighbour geometry favoured by 
the quadrupole-quadrupole interactions encourages the 
formation of two-dimensional square networks. However, 
this tendency is frustrated by the LJ part of the poten- 
tial, which drives the structure towards compact, highly- 
coordinated arrangements with polytetrahedral or icosa- 
hedral packing [25[. This competition produces a series 
of structural motifs that partially satisfy the two oppos- 
ing trends. Fig. [3] summarises the structural maps of 
SQat and LQn as a function of the number N of par- 
ticles and the strength Q* of the quadrupole moment. 
Some structures are difficult to classify in an unambigu- 
ous or meaningful way, and such combinations of N and 
Q* have been left blank in the figure for clarity. 

For a sufficiently weak quadrupole moment, the 
global minimum must be close to the LJ global mini- 
mum structure, but with slight distortions induced by 
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FIG. 3: Structural maps for (a) the square and (b) the linear 
quadrupole: relaxed LJ structure (filled circle), sheet (filled 
square), stacked triangular antiprisms (open upward trian- 
gle), decorated triangular antiprisms (filled upward triangle), 
stacked square antiprisms (open square), filled stacked pen- 
tagonal antiprisms (plus), hollow shell (open circle), filled 
shell (filled downward triangle), lattice- like (cross), decahe- 
dral core (filled diamonds), star (star). 




FIG. 4: Structures discussed in the text, (a) SQi6 sheet at 
Q* = 2, (b) LQ16 sheet at Q* = 1.5, (c) SQ13 cuboctahedron 
at Q* = 0.4, (d) SQ 24 filled shell at Q* = 0.6, (e) SQ 24 
decorated triangular antiprisms at Q* = 0.4, (f) LQ12 star at 
Q* = 1.025. 



the quadrupoles. However, effectively confining the 
quadrupoles to an icosahedral framework for SQ13 causes 
the quadrupoles to experience severe frustration, akin 
to that found in geometrically frustrated magnets [3l| . 
Hence, a given arrangement of the particles can corre- 
spond to multiple potential energy minima in the orien- 
tational part of configuration space, introducing a new 
source of complexity to the PES. The number of such 
isomers of the LJ structure generally increases with Q*, 
but also depends sensitively and non-monotonically on 
N. For the near-icosahedral SQ13 there are two distinct 
isomers at Q* = 0.025, while for SQ19, where the global 
minimum is based on two interpenetrating icosahedra for 
this value of Q*, we located 23 distinct arrangements 
by quenches of the LJ structure starting from random 
quadrupole orientations. This figure rises to several hun- 
dred for SQ19 at Q* = 0.1. Although these searches 
are not definitively exhaustive, the rapid increase in the 
number of isomers illustrates the roughness of the PES 
with respect to the orientational coordinates. 

In the opposite limit of large Q*, a two-dimensional 
sheet always emerges. The sheets for the square num- 
ber N = 16 are shown in Fig. [H The twist seen in 
LQ4 is continued as the sheet grows, while the sheets of 
square quadrupoles are always planar. The sheet grows 
by adding particles at adjacent sites of the extended 
square lattice, with the dimensions of the lattice adapt- 
ing to maximise the number of T-shaped pairs in the 
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FIG. 5: Part of some growth sequences for the square 
quadrupole. (a) stacked triangular antiprisms, (b) stacked 
square antiprisms, (c) hollow shells. 



first instance. It is often possible to achieve the maxi- 
mum number of such pairs in more than one way on a 
square lattice, and next-nearest neighbour interactions 
then come into play. Hence, there can be close compe- 
tition between structures even in the strong quadrupole 
limit. In contrast, the Stockmayer potential always has 
an unambiguous optimal structure consisting of a planar 
ring of head-to-tail dipoles in the strong dipole limit [lo[ . 
The lower boundary of the quadrupole sheet on the struc- 
tural map moves (non-monotonically) to higher Q* as TV 
increases, because a larger number of LJ pair interactions 
must be disrupted to create the sheet. The widening re- 
gion between the relaxed LJ cluster and the sheet is oc- 
cupied by structures that strike a compromise between 
high-coordination number and sheet-like arrangements. 



C. The 13-particle cluster 

The 13-particle cluster, which for the pure LJ poten- 
tial has a global minimum consisting of a centred icosa- 
hedron with point group J^, provides a good illustration 
of the sequence of changes driven by the quadrupole in 
small clusters. For SQ13 with small Q*, the quadrupole 
of the central particle aligns itself perpendicular to one 
of the icosahedral C2 axes. The resulting small distor- 
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FIG. 6: The 13-particle icosahedron (point group Ih) viewed 
along a C3 axis and the two LQ13 isomers at Q* = 0.025 of 
lowest energy. Both belong to point group #6. 



tions of the particle positions lower the symmetry to 
C<ih- The quadrupole-quadrupole interactions are highly 
frustrated when confined to the vertices of the distorted 
icosahedron and at Q* = 0.35, the global minimum 
switches to a slightly distorted centred cuboctahedron 
(Pid) — a fragment of face-centred cubic lattice (Fig. H^). 
This structure maintains the 12-fold coordination of the 
central particle but the shell consists of squares and tri- 
angles, with quadrupole arrangements closer to those in 
the right-hand panels of Fig.[2j At Q* = 0.525, a second 
change occurs, to a stack of face-sharing triangular an- 
tiprisms (illustrated in the second panel of Fig. [5^), which 
can also be regarded as a narrow tube if viewed down 
the three-fold axis. This structure sacrifices the high co- 
ordination of a more spherical lattice-like fragment for 
an elongated arrangement, in which the quadrupoles are 
better aligned. The switch to the sheet structure then 
takes place in two steps. First, at Q* = 0.8 a 3 x 4 sheet 
arises with the thirteenth particle in the same plane, 
bridging the central bond of a long side. At Q* = 1.15, 
the triangular face becomes too unfavourable, and the 
sheet adopts a 4 x 4 square with three of the corners 
missing. 

A related sequence emerges with Q* for the cluster of 
13 linear quadrupoles. Like the dipole moment in the 13- 
particle Stockmayer cluster, the axis of the quadrupole on 
the central particle in LQ13 selects one of the C3 axes of 
the icosahedron, here reducing the symmetry from 1^ to 
Sq. However, there are four other icosahedral minima dif- 
fering by the quadrupole orientations. The two isomers 
with the lowest energy are illustrated in Fig. [6j With in- 
creasing Q*, the cluster eventually passes to the stacked 
triangular ant ipr isms, but does so through a rather amor- 
phous structure, unlike the highly symmetric cuboctahe- 
dron seen in SQ13. The linear quadrupole can tolerate a 
considerably higher Q* before switching to the sheet than 
can the square quadrupole — an observation that holds for 
all N studied here. 



D. Structural families and growth sequences 

A number of structural families emerge in LJQ clusters 
and persist over some range of Q* and N. The stacked 
triangular antiprisms described above for 13 particles are 
seen for SQat with 9 < N < 15. When N is not a multi- 



ple of three, the end particles form an incomplete layer, 
giving rise to a simple growth sequence, part of which 
is depicted in Fig. [5^. This family does not continue 
indefinitely as the global minimum, but is replaced for 
16 < N < 21 by stacked square antiprisms (Fig. (5}d). 
The diamond-like faces on the surface of the square-based 
structure are flatter and closer to the ideal tetramer ar- 
rangement than those on the surface of the triangular 
antiprismatic stack. 

Increasing the circumference of these stacks by another 
particle to make pentagons makes the diameter of the 
structure large enough to accommodate a line of parti- 
cles down the centre of the stack, giving a filled tube-like 
arrangement. The additional contacts provided by the 
central line make such structures competitive at lower 
Q* than the stacked squares. In fact, pentagonal stacks 
arise in two forms in the structural map. The diamond 
symbols in Fig. [3](a) indicate clusters built around a dec- 
ahedron, which contains a pentagonal prism. For pair 
potentials with a preferred nearest-neighbour separation, 
decahedral structures are less strained than icosahedral 
ones and are seen in the global minima of short-ranged 
isotropic potentials [32[ and experimentally in metal clus- 
ters j33| . In LJQ clusters, the square faces are favoured 
by the quadrupolar interactions, and the structure can 
grow by building additional partial layers around the 
five-fold axis. Twisting the pentagonal layers gives the 
second type of pentagonal stack structure, filled pentag- 
onal antiprisms, which appear on the map for SQat with 
N > 20. These tube-like stacks can be regarded as two- 
dimensional sheets in which two opposite edges have been 
joined, thereby exchanging the energetic cost of an ex- 
posed edge for the penalty of curving the sheet. This 
trade-off is analogous to the formation of closed rings of 
dipoles [20[. Larger tubes have been observed to assem- 
ble spontaneously in the finite-temperature simulations 
of the LJQ fluid by Van Workum and Douglas [211 ]. 

Tubes can dispose of their remaining exposed edges by 
also closing the ends to make a shell. We observe hollow 
shells over a range of Q* in SQ22 and larger. For this class 
of structures, certain values of N give rise to a structure 
of high symmetry. The first of these is SQ24, illustrated 
in the central panel of Fig. [5b, which has perfect Oh octa- 
hedral symmetry. The shell normally grows by insertion 
of a particle into the surface, causing a distortion of the 
ideal triangular and square faces, but occasionally by the 
addition of an edge-bridging particle, as shown in Fig.[5b. 
The shell seems to be a permanent feature of larger SQat 
clusters. We have followed it as far as TV = 36, which 
forms an elongated shell of Dsd symmetry with triangular 
faces at the ends and antiprismatically stacked hexagons 
along the body. 

If Q* is not sufficiently large, the shell is energetically 
penalised for its shortage of LJ nearest neighbour pairs. 
However, a large number of such pairs can be obtained 
by placing a few particles inside the shell. Hence, the 
hollow shell is typically preceded by a filled shell or filled 
pentagonal tube in the structural map, Fig. [3fci. A shell 
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encapsulating two particles is shown in Fig. [4]1 for SQ24; 
compare the hollow shell for this cluster in the central 
panel of Fig. Et. 

Similar energetic compromises produce mixtures of 
structures that have already been described. For exam- 
ple, the transition from the LJ structure to the filled 
tubes and shells is sometimes bridged by decorated ver- 
sions of the stacked triangular antiprisms, where the 
stack has been surrounded by a new layer. This arrange- 
ment is illustrated in a view down the three-fold axis for 
SQ24 in Fig. [4^. The characteristic network of square 
faces for the shell is beginning to emerge on the outside 
of these structures, but they maintain a larger number of 
LJ pairs than the filled shell. 

The linear quadrupole tends not to give rise to hollow 
global minima. Although tube-like and shell-like struc- 
tures do appear, they are collapsed into what would be 
the central space in the square quadrupole equivalents. 
This ability to distort, or inability to support a hollow 
interior, is a result of the axial symmetry of the linear 
quadrupole. A square array of linear quadrupoles, such 
as the one depicted in Fig. [3t>, can fold and twist along 
one of its diagonals without severely disrupting the T- 
shaped nearest-neighbour interactions either side of the 
fold. This is not true of an array of square quadrupoles, 
such as that in Fig. [4^, in which each quadrupole defines 
a plane and not just a line. By collapsing inwards, the 
clusters of linear quadrupoles gain favourable interactions 
between opposite sides of the structure that would other- 
wise be held far apart. However, the collapsed structures 
are often rather amorphous, making them hard to clas- 
sify or describe in a helpful way. For this reason, the 
structural map in Fig. [3]d extends only to LQ20. 

A distinctive and reproducible feature of the linear 
quadrupole clusters is a family of structures with a star- 
like organisation of the particles and a gear-wheel ar- 
ray of quadrupole axes, which are scattered around the 
structural map (Fig. [3Jd). Again, for particular values of 
iV, the cluster can achieve a high symmetry that may 
be based on a three-fold or four-fold principal symmetry 
axis. An example belonging to point group Sq is shown 
in Fig. [4f. However, the stability of these morphologies 
is strongly correlated with the number of particles; the 
stars are not observed away from the values of N that 
allow the symmetry to be completed. 



IV. CONCLUDING REMARKS 

The survey of putative global minima presented here 
shows that the competition between isotropic attractive 
forces and highly directional quadrupole-quadrupole in- 
teractions gives rise to a wide variety of structural motifs. 
These include elongated tube-like structures, hollow and 
filled shells, stars, and extended sheets. Some unusual 
point groups are represented in this collection. 

Global optimisation is most challenging when the com- 
peting influences are closely balanced, i.e., for interme- 



diate strengths of the quadrupole in this work. Ob- 
taining reproducible lowest-energy structures was signif- 
icantly more difficult for the quadrupolar potential than 
for equivalent dipolar Stockmayer potential [20] . This ob- 
servation, together with the multiplicity of minima that 
were found to have virtually identical positions but dif- 
ferent orientations of the quadrupoles, hints at a complex 
potential energy surface in certain parts of the (N,Q*) 
parameter space. Confirmation and further exploration 
of this complexity would require a more comprehensive 
analysis of the energy landscape [34| . The landscape ap- 
proach would provide information not only on the num- 
ber of competing structures, but also on the barriers sep- 
arating them and the rearrangement mechanisms that 
interconvert them. This information would provide a 
starting point for investigating the thermal stability and 
dynamic properties of the clusters, as well as the routes 
by which they might self- assemble. 

The various families of structures have their own 
"magic" numbers at which a particular shape is com- 
plete and the landscape is probably minimally frustrated 
jlil |35| ]. Such numbers are well known for a variety 
of simpler interatomic potentials, and often correspond 
to the completion of successive icosahedral shells [36| 
at N = 13, 55, 147.... The special stability associ- 
ated with these sizes, combined with kinetic accessibility 
[l4|, can lead to prominent features such as experimen- 
tal abundance [37[. In the present work, the tube-like 
structures of stacked triangular and square antiprisms 
achieve completed layers for multiples of three and four 
particles, respectively, while a shell can be elongated by 
the insertion of a complete hexagonal ring. In the strong 
quadrupole regime, sheets adopt defect-free squares when 
N is a square number. It would be interesting to inves- 
tigate whether these "perfect" structures are especially 
stable and self-assemble efficiently, as for magic number 
Lennard- Jones clusters [38]. 

We have seen that quadrupole-quadrupole interac- 
tions favour the formation of extended two-dimensional 
structures with four-fold coordination of the particles. 
In contrast, dipole-dipole interactions lead to extended 
pseudo one-dimensional chains, while isotropic attraction 
drives the structure towards compact three-dimensional 
arrangements. By careful balancing of the multipolar in- 
teractions, it should therefore be possible to exert consid- 
erable control over the structures that self-assemble out 
of multipolar particles with isotropic core interactions. 

Briefly considering bulk phases, rather than finite clus- 
ters, such control could be useful in adjusting the net- 
working properties of colloidal gels. For example, it has 
recently been shown that dipolar colloids can be encour- 
aged to form more interconnected gel-like networks by a 
slight extension of the dipole [39[ . From studies of models 
of patchy spheres with fixed maximum valency it is now 
known that the average coordination number of the par- 
ticles in a gel has important consequences for the struc- 
ture of the gel and for the underlying phase behaviour 
of the fluid from which it forms [40[. The present work 
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suggests that the average coordination number could be 
finely tuned either by adding a weak point quadrupole 
to point dipolar particles, or by using a mixture of dipo- 
lar and quadrupolar spheres. Hence, multipolar particles 
could be an appealing alternative to patchy colloids for 
realizing and exploring reversible gels [411 ]. 
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